Reactor On- Off Antineutrino Measurement with KamLAND 
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The recent long-term shutdown of Japanese nuclear reactors has resulted in a significantly reduced reactor 
flux at KamLAND. This running condition provides a unique opportunity to confirm and constrain back- 
grounds for the reactor 77^ oscillation analysis. The data set also has improved sensitivity for other Ve signals, 
in particular 77e's produced in /3-decays from ^•^^U and -^^-^Th within the Earth's interior, whose energy spec- 
trum overlaps with that of reactor 77e's. Including constraints on 613 from accelerator and short-baseline reactor 
neutrino experiments, a combined three-flavor analysis of solar and KamLAND data gives fit values for the os- 



cillation parameters of tan^ 6*12 = 0.436lo;o25' ^^21 = 
Assuming a chondritic ThAJ mass ratio, we obtain llGj 
-+0.8 



' eV^ and sin^ 6, 



0.0231 



e events from ^-^^U and '^^'^Th, corresponding to 
a geo Ve flux of 3.4j;g g x 10^ cm~^s~'^ at the KamLAND location. We evaluate various bulk silicate Earth 
composition models using the observed geo Ve rate. 

PACS numbers: 14.60.Pq, 28.50.Hw, 91.35.-x, 91.67.Qr 
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I. INTRODUCTION 

The Kamioka Liquid-scintillator Antineutrino Detector 
(KamLAND) demonstrated the oscillatory nature of neutrino 
flavor transformation by observing electron antineutrinos (Fg) 
with energies of a few MeV from nuclear reactors typically 
180 km away Ol]. Following the Fukushima nuclear accident 
in March 2011, the entire Japanese nuclear reactor industry, 
which generates >97% of the reactor Ve flux at KamLAND, 
has been subjected to a protracted shutdown due to a review 
of nuclear safety standards. This unexpected situation allows 
for a reactor on-off study of backgrounds for the KamLAND 
neutrino oscillation analysis. 

The reactor-off data also yield improved sensitivity for 
Fe's produced by other sources. Previously, we used the 
KamLAND data to search for geoneutrinos, 77e's produced 
in /3-decays from primordial radioactivity within the Earth's 
interior The ^asjj ^jj^j (jgc^y chains emit Fg's with en- 

ergies below 3.4 MeV, so reactor events with similar ener- 
gies pose a background for this signal. Despite having a high 
reactor T7g background, KamLAND performed the first exper- 



imental study of geo Fe's from the decay chains of 238 y ^jj^j 
^■^^Th Later the geo signal at KamLAND was used 
to estimate our planet's radiogenic heat production and con- 
strain composition models of the bulk silicate Earth (BSE, the 
Earth's region outside its metallic core). In particular it was 
found that fully -radiogenic Earth models are disfavored [jst]. 
The Borexino experiment at Gran Sasso also reported a posi- 
tive observation of geo Fe's jH]. 

In this article, we present improved reactor neutrino oscilla- 
tion results and geo F,. flux measurements including the recent 
reactor-off period. For the reactor rate estimate, we also ap- 
ply new evaluations of reactor antineutrino emission spectra, 
as well as constraints on oscillation parameters from acceler- 
ator and short-baseline reactor neutrino oscillation measure- 
ments. 



II. NEUTRINO OSCILLATION 

Neutrino oscillation is well established by experimental 
studies of solar, reactor, atmospheric, and accelerator neutri- 
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nos. KamLAND observes Ve's from many reactors at a flux- 
weighted average distance of 180 km, providing optimal sen- 
sitivity for the LMA-MSW i^i-i^2 mixing solution to the solar 
neutrino problem. For the length scale relevant to reactor 
oscillation at KamLAND, the three-flavor survival probability 
(P^^), including matter effects, may be approximated as 



cos^ Oi.P^: + sin^ 
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The two neutrino survival probability P'^^^ has the same form 
as the survival probability in matter for vi-V2 mixing but with 



the electron density {N^) modified: 
is given by 



iVeCOs2 0i3 H. It 
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where L is the distance from the source to the detector, E^, is 
the energy, and 9i2m and Am^u^.j are the matter-modified 
mixing angle and mass splitting defined by 



sin^ 26 



sin^ 26*12 
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(cos 26*12 - A/Am|i)2 + sin^ 26*12 



,(3) 



^^2iM = Am^i^(cos26'i2 - yl/Am|i)2 + sin^ 26*12 -(4) 

The parameter A = ±2\/2GFNeEi, has a negative sign for 
antineutrinos; Gf is the Fermi coupling constant. 

Recently, accelerator and short-baseline (~lkm) reactor 
experiments have demonstrated that£i 3 is non-zero, and have 
measured it with high precision ll6l-[l0ll. An analysis incorpo- 
rating this new ^^13 constraint will improve the determination 
of the other oscillation parameters. 



III. GEONEUTRINO FLUX AT KAMLAND 

While the mechanical properties of the Earth's interior are 
well established, its composition, including its radiochemical 
content, remains uncertain. Decays of uranium (U), thorium 
(Th), potassium (K) and their progeny generate heat. Depend- 
ing on their abundance and distribution within the Earth, these 
decays may be an essential heat source for generating Earth 
dynamics. A leading BSE model ifllll based on measured el- 
emental abundances of chondritic meteorites and mantle peri- 
dotites predicts a radiogenic heat production of 8 TW from 
the 238u decay chain, 8TW from the ^"^^Th decay chain, and 
4 TW from fT2^ . This would account for nearly half of the 
heat dissipation rate from the Earth's surface, which a recent 
analysis finds to be 47 ± 2 TW fill. 

The energy spectrum of '^"K neutrinos falls entirely below 
the 1.8 MeV energy threshold for the inverse /3-decay reaction 
by which KamLAND observes antineutrinos, rendering these 
decays invisible to KamLAND. However, the ^^^U and ^^^Th 
decay chain Pg's extend above this threshold with distinct en- 
ergy distributions, making possible a direct measurement of 
the individual ^38^ ^jj^j 232'pjj contributions. 
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FIG. 1 : Schematic diagram of the KamLAND detector. The shaded 
region in the liquid scintillator indicates the volume for the Ve anal- 
ysis after the inner balloon was installed. 



The geo Ve flux at the KamLAND detector can be calcu- 
lated from the isotope abundances ai{r') for each isotope i at 
source positions r' by integrating over the entire Earth, 



dE, 



dE,, 
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where r is the detector position, Ai is the decay rate per unit 
mass, dui (E^)/ dEi, is the Ve energy spectrum for each mode 
of decay, a; (r') is the isotope mass per unit rock mass, p(r') is 
the rock density, and Pee{E^, \ f~ r'|) is the Ve survival prob- 
ability given by Eq. ([T]) with L — \ f~r'\. Given the measured 
values of neutrino oscillation parameters and the energy range 
of detectable geo Fg's, the integration over the volume of the 
Earth averages over the second sine function in Eq. (|2}, allow- 
ing the approximation, 



cos (^13 



1 - i sin^ 2012 



sm 6'i3. 
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In Eq. ([6]) we have neglected matter effects, which modify 
the survival probability by <1% lfl4ll . From a global analysis 
of neutrino oscillation data involving solar, accelerator, and 
reactor neutrinos, including the present KamLAND data, we 
obtain Pee — 0.551 ± 0.015. The less-than 3% uncertainty 
in Pee is negligible compared to the statistical uncertainty of 
KamLAND 's current geo Ve flux measurement. 



IV. THE KAMLAND EXPERIMENT 

KamLAND is located in Gifu Prefecture, Japan, under 
Mount Ikenoyama at a depth of ~'2700m water-equivalent. 
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The primary volume consists of 1 kton of ultra-pure liquid 
scintillator (LS), which comprises the neutrino interaction tar- 
get (Fig. [TJ. The LS is contained in a 13-m-diameter spher- 
ical balloon made of 135-/^jm-thick transparent nylon/EVOH 
(ethylene vinyl alcohol copolymer) composite film. The bal- 
loon is suspended in non-scintillating purified mineral oil con- 
tained inside an 18-m-diameter stainless steel tank. The LS 
consists of 80% dodecane and 20% pseudocumene (1,2,4- 
trimethylbenzene) by volume, and 1.36 ± 0.03 g/liter PPO 
(2,5-diphenyloxazole) as a fluor The scintillation light is 
viewed by an array of 1325 specially-developed fast 20-inch- 
diameter photomultiplier tubes (PMTs) masked to 17-inch 
diameter, and 554 older unmasked 20-inch PMTs, provid- 
ing 34% solid-angle coverage in total. This inner detector 
(ID) is surrounded by a 3.2 kton water-Cherenkov outer de- 
tector (OD) that serves as a cosmic -ray muon veto counter 

In September 2011, the KamLAND-Zen neutrinoless dou- 
ble beta-decay search was launched ifTsll . This search makes 
use of KamLAND's extremely low background and suspends 
a /?/? source, 13 tons of Xe-loaded liquid scintillator (Xe-LS), 
in a 3.08-m-diameter inner balloon (IB) at the center of the 
detector, as shown in Fig. [1] To avoid backgrounds from the 
IB and its support material, the analysis reported here is 
restricted to events occurring well outside the IB. 

Electron antineutrinos are detected through the inverse (3- 
decay reaction, + p ^ e+ + n, which yields a delayed co- 
incidence (DC) event pair signature that provides a power- 
ful tool to suppress backgrounds. The prompt scintillation 
light from the gives a measure of the incident Ve energy, 
c:^ Ep + En + 0.8 MeV, where Ep is the sum of the e+ 
kinetic energy and annihilation 7 energies, and En is the av- 
erage neutron recoil energy, 0(10 kcV). The mean time for 
capture of the neutron in the LS is 207.5 ± 2.8 fis lfl6ll . The 
scintillation light from the capture 7 constitutes the delayed 
event of the DC pair 



V. ANTINEUTRINO CANDIDATE EVENT SELECTION 

The data reported here are based on a total live-time of 
2991 days, collected between March 9, 2002 and Novem- 
ber 20, 2012. The data set is divided into three periods. 
Period 1 (1486 days live-time) refers to data taken up to May 
2007, at which time we embarked on a LS purification cam- 
paign that continued into 2009. Period 2 (1154 days live- 
time) refers to data taken during and after the LS purifica- 
tion campaign, and Period 3 (35 1 days live-time) denotes the 
data taken after installing the IB. We removed periods of low 
data quality and high dead time that occurred during LS pu- 
rification and KamLAND-Zen IB installation. The LS pu- 
rification reduced the dominant Period 1 background for P^'s, 
^'^0(0;, n)^^0 decays, by a factor of ^20. The high-quality 
data taken after LS purification accounts for 50% of the to- 
tal live-time. Using a spherical fiducial scintillator volume 
with 6.0 m radius, the number of target protons is estimated 
to be (5.98 ± 0.13) x lO^^, resulting in a total exposure of 
(4.90 ± 0.10) X lO'^^ target-proton-years. The reduced fidu- 
cial volume in Period 3 is accounted for in the detection ef- 



ficiency; it contributes negligible additional fiducial volume 
uncertainty for Period 3. 

Event vertex and energy reconstruction is based on the 
timing and charge distributions of scintillation photons 
recorded by the ID PMTs. The reconstruction is caU- 
brated with ^"Co, ^^Ge, "^"^Ug, ^^Zn, ^^lAm^Be, "7cs, and 
210pq13(^ radioacti ve sources . The achieved vertex resolu- 
tion is ^I2cm/ ^yE{McV), and the energy resolution is 
6.4% / y/ E{MeV). The nonlinear, particle-dependent con- 
version between deposited (real) energy and KamLAND's 
prompt energy scale is performed with a model incorporat- 
ing Birks quenching and Cherenkov emission. The model 
parameters are constrained with calibration data, and con- 
tribute a 1.8% systematic uncertainty to the measured value 
of Ato2i- Using calibration data taken throughout the fiducial 
volume during Period 1, we find that the deviation of recon- 
structed vertices from the actual deployment locations is less 
than 3 cm. Incorporating a study of muon-induced ^^B/^^N 
decays ifisll . the fiducial volume uncertainties are 1.8% for the 
pre-purification data and 2.5% for the post-purification data. 

For the DC event pair selection, we apply the following 
series of cuts: (i) prompt energy: 0.9 < £'p(MeV) < 8.5; 

(ii) delayed energy: 1.8 < Ed{McV) < 2.6 (capture 
on p), or 4.4 < iJd(MeV) < 5.6 (capture on ^^C); 

(iii) spatial correlation of prompt and delayed events: 
Ai?<2.0m; (iv) time separation between prompt and 
delayed events: 0.5 < AT{fis) < 1000; (v) fiducial volume 
radii: Rp, < 6.0 m; (vi) and for Period 3, delayed vertex 
position: i?d > 2.5m and pd > 2.5m, Zd > Om (vertical 
central cylinder cut at the upper hemisphere) to eliminate 
backgrounds from the KamLAND-Zen material. To maxi- 
mize the sensitivity to Fe signals, we perform an additional 
event selection designed to suppress accidental coincidence 
backgrounds from radioactive contaminants in the detector 
while maintaining high efficiency for Fe's. First, using a 
combination of Monte Carlo, data-driven, and analytical 
methods, we constructed a probability density function (PDF) 
for the I7g signal (fj^J and accidental (face) coincidence 
events. The PDF is based on the six cut parameters {Ep, E^, 
AR, AT, Rp, i?d)- For each candidate pair, we calculate the 
discriminant L ~ J^^+j — determine a selection value, 
Lcut{Ep), to maximize the figure-of-merit for 
prompt energy intervals of 0. 1 MeV. In the figure-of-merit, 
S is the number of the expected signal events assuming an 
oscillation-free reactor spectrum and the geo Ve fluxes pre- 
dicted by UtIi . Bacc corresponds to the number of accidental 
background events, as measured using an out-of-time delayed 
coincidence window selection (10 ms < AT < 20 s). The 
selection efficiency is calculated via Monte Carlo from the 
ratio of selected Ue's to the total number of generated Fg's in 
i? < 6 m. The systematic uncertainty is evaluated using ^^Ge 
and ^^^Am^Be source calibrations as discussed in [OJ. The 
total number of events passing all selection criteria is 261 1. 

The reactor fluxes can be predicted from reactor operation 
records, which are provided to the KamLAND Collaboration 
by a consortium of Japanese electric power companies, and 
include the thermal power variation as well as fuel exchange 
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TABLE I: Estimated backgrounds for Ve in tlie energy range between 0.9 MeV and 8.5 MeV after event selection cuts. 



Background 


Period 1 


Period 2 


Period 3 


All Periods 






(1486 days) 


(1154 days) 


(351 days) 


(2991 days) 


1 


Accidental 


76.1 ± 0.1 


44.7 ±0.1 


4.7 ±0.1 


125.5 ± 0.1 


2 


^Li/^He 


17.9 ± 1.4 


11. 2± 1.1 


2.5 ±0.5 


31.6 ± 1.9 




^^C(a, n)^'^Og.s., elastic scattering 


160.4 ± 16.4 


16.5 ±3.8 


2.3 ± 1.0 


179.0 ± 21.1 




. i3C(a, n)i«Og.s., i2C(n, n')'^C* (4.4 MeV 7) 


6.9 ± 0.7 


0.7 ± 0.2 


0.10 ±0.04 


7.7 ± 0.9 




f ^^C(a, 71)^*^0*, 1st e.s. (6.05 MeV e+e") 


14.6 ± 2.9 


1.7 ±0.5 


0.21 ± 0.09 


16.5 ± 3.5 


*l 


"C(q, n)^'^0*, 2nd e.s. (6.13 MeV 7) 


3.4 ± 0.7 


0.4 ±0.1 


0.05 ± 0.02 


3.9 ± 0.8 


5 


Fast neutron and atmospheric neutrino 


< 7.7 


<5.9 


< 1.7 


< 15.3 


Total 


279.2 ±22.1 


75.2 ±7.6 


9.9 ±2.1 


364.1 ±30.5 




FIG. 2: Time evolution of expected and observed rates at KamLAND for Fe's with energies between (a) 0.9 MeV and 2.6 MeV and (b) 
2.6 MeV and 8.5 MeV. The points indicate the measured rates in a coarse time binning, while the curves show the expected rate variation for 
reactor Fe's (black line), reactor F^'s ± backgrounds (colored line), and reactor Fe's ± backgrounds + geo Fe's (gray line). The geo rates 
are calculated from the reference model iUtIi . The vertical bands conespond to data periods not used in the analysis. In the right panel of (a), 
the data are grouped according to periods of similar expected reactor + background rates, as denoted by the colored bands. The observed 
event rate for each group is plotted at the exposure-weighted expected event rate within the group. The efficiency-corrected best-fit value of 
the geo Vc rate from the full spectral analysis (dashed line), its la error (shaded region), and the model expectation (gray line) are drawn for 
comparison. The contribution of geo Fe's in (b) is negligible. The oscillation parameters used to calculate the expected reactor Fe rate are the 
best-fit values from the global oscillation analysis: tan^ 6I12 = 0.436tJ^:Jj2^, Amii = 7.53lE];i^ x 10~^ eV^ and sin^ 6I13 = 0.023l!5;^^2- 



and reshuffling data for afl Japanese commercial reactors. The 
thermal power generation used for the normalization of the 
fission rates is measured to within 2%. Only four isotopes 
contribute significantly to the emission spectra; the relative 
fission yields, averaged over the entire live-time period for this 
result, are (0.567 : 0.078 : 0.298 : 0.057) for (^a^u : ^ssy : 
239p^ . 24ip|^-j^ respectively. A recent recalculation of the 77e 
spectra per fission of these isotopes introduced a ^3% upward 



shift lfl9[|20ll relative to the previous standard calculation ll2ll 
I22I1 . causing past measurements at short-baselines to appear 
to have seen fewer Pe 's than expected. It has been speculated 
that this so-called Reactor Antineutrino Anomaly may be due 
to some systematic uncertainty or bias, or could potentially 
be due to oscillation into a heavy sterile neutrino state with 
Am ~ leV^ iH. To make our analysis insensitive to these 
effects, the normahzation of the cross section per fission for 
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each reactor is adjusted to reproduce the Bugey4 result 112411 : 



\^lreac. ^ BiigeyA 



Bugey4: 



)(^). (7) 



where ai is the fractional fission rate of the isotope i. The 
contribution from Korean reactors, based on reported elec- 
tric power generation, is estimated to be (4.9 ± 0.5)%. The 
contribution from Japanese research reactors and all other re- 
actors around the world is (1.1 ± 0.6)%. The levels of the 
long-lived, out-of-equilibrium fission products ^"Sr, ^°^Ru, 
and "-^Ce IH are evaluated from the history of fission rates 
for each isotope and are found to contribute an additional 
(0.7 ± 0.3)%. Applying the selection cut efficiency, we ex- 
pect a total of 3564 ± 145 events from reactors in the absence 
ofVe disappearance. 

A calculation of the geo Ve flux at KamLAND based on 
the reference Earth model of iUtI gives an expected 109 and 
27 geo Ve events from U and Th, respectively. Since the es- 
timation of the geo yield is highly model-dependent, the 
event rates from the U and Th decay chains are not constrained 
in the oscillation analysis. Only the prompt energy spectral 
shapes, which are independent of the Earth model, are used to 
constrain their contributions. A possible contribution from a 
hypothetical reactor- source at the Earth's center, motivated 
by ll26ll and investigated in |40 and jsfli is neglected as a back- 
ground in the fit for the oscillation parameters and geoneutrino 
fluxes, but is addressed briefly below as an independent signal. 

In Period 1, the dominant background is the ^^C(a, n)^^0 
reaction, generated from the a-decay of ^^"Po in the LS. 
The neutrons in this reaction are produced with energies 
up to 7.3 Me V, but the visible energy is quenched to be- 
low 2.7 Me V. Accounting for the energy-dependent effi- 
ciency of the Lcut{Ep) selection, the estimated number of 
"C(a, n)^^0 background events is 207.1 ± 26.3 in the en- 
ergy region 0.9 < i?p(McV) < 8.5. The accidental back- 
ground, which dominates in Periods 2 and 3, is measured with 
an out-of-time delayed coincidence window from 10 ms to 
20 s to be 125.5 ±0.1 events. Including smaller contributions 
from cosmogenically produced radioactive isotopes, fast neu- 
trons produced by cosmic-ray muons, and atmospheric neu- 
trinos, the total background is estimated to be 364.1 ± 30.5 
events. The backgrounds are detailed in Table J] 
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FIG. 3: Prompt energy spectrum of Ve candidate events above the 
0.9 MeV energy threshold (vertical dashed line) for each data tak- 
ing period. The background, reactor and geo Ve contributions are 
the best-fit values from a KamLAND-only analysis. The prompt en- 
ergy spectra of Ve candidate events in the low-energy region are also 
shown in the inset panels with a finer binning. The top panel shows 
the energy-dependent selection efficiency curves for each period. 



VI. ANTINEUTRINO MEASUREMENT AND 
OSCILLATION ANALYSIS 

To extract the neutrino oscillation parameters and geoneu- 
trino fluxes, 77e candidates are analyzed with an unbinned 
maximum-likelihood method incorporating the event rate and 
the prompt energy spectrum shape, including their time vari- 
ation, in the range 0.9 < £'p(McV) < 8.5. The is defined 
by 

= X?atc(^'l2,^'l3, ATO21,^^BGl-^5,^^u!Th'"l^4) 

-2 In ishapo (6*12 , Oi3 , Am2i , A'^BGl^.S , ^S'^Th 7 ai->4) 
+ xIg(^BG1->5) + Xsyst("l->4) 

+xLi(^i2,^i3,A™2j. (8) 



The terms are, in order: the contribution for (i) the time- 
varying event rate, (ii) the time-varying prompt energy spec- 
trum shape, (iii) a penalty term for backgrounds, (iv) a penalty 
term for systematic uncertainties, and (v) a penalty term for 
the oscillation parameters. N^°^y^ are the flux normalization 
parameters for U and Th geo X^e's, and allow for an Earth- 
model-independent analysis. iVeci^s are the expected num- 
ber of backgrounds, and are allowed to vary in the fit but are 
constrained with the penalty term (iii) using the estimates de- 
scribed in the preceding section and listed, with the corre- 
sponding index, in Table U ai^4 parametrize the uncertain- 
ties on the reactor spectrum, the energy scale, the event 
rate, and the energy-dependent detection efficiency; these pa- 
rameters are allowed to vary in the analysis but are constrained 
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FIG. 4: Allowed regions projected in the (tan 612, Am^i) plane, 
for solar and KamLAND data from the three-flavor oscillation anal- 
ysis for (a) 6*13 free and (b) 613 constrained by accelerator and short- 
baseline reactor neutrino experiments. The shaded regions are from 
the combined analysis of the solar and KamLAND data. The side 
panels show the Ax'^-profiles projected onto the tan^ 612 and Am2i 
axes. 



by term (iv). Table Ull summarizes the systematic uncertainties 
on Am2i and the expected event rate of reactor Fe's. The 
overall rate uncertainties for Period 1 and for Periods 2 and 3 
are 3.5% and 4.0%, respectively. Systematic uncertainties 
are conservatively treated as being fully correlated across all 
data taking periods. The penalty term (v) optionally provides 
a constraint on the neutiino oscillation parameters from so- 
lar ||2%|3l|l, accelerator (T2K t^, MINOS lltl), and short- 
baseline reactor neutrino experiments (Double Chooz 
Daya Bay §], RENO fioll ). 
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FIG. 5: Ratio of the observed Ve spectrum to the expectation for 
no-oscillation versus Lo/E for the KamLAND data. Lo — 180 km 
is the flux-weighted average reactor baseline. The S-u histogram is 
the best-fit survival probability curve from the three-flavor unbinned 
maximum-likelihood analysis using only the KamLAND data. 



Figure|2]plots the time variation for the rates of reactor T^e's, 
geo Fg's, and backgrounds for the three data taking periods, 
assuming the best-fit oscillati on p arameters, and geo Ve fluxes 
from the reference model of iflTII . Also drawn are the correla- 
tions between the measured and expected best-fit event rates, 
which should fit to a line with unit slope and zero offset in the 
absence of geo Fe's. The vertical displacement of the trend 
for events below 2.6 MeV is attributed to the contribution of 
geo 77e's. 

Figure [3] shows the prompt energy spectra of Ve candidate 
events for each period. The reduction of the ^^C(a, n)^^0 
background in Period 2 and of reactor Fe's in Period 3 can 
clearly be seen. For the three-flavor KamLAND-only anal- 
ysis (Xosci = 0), the fit oscillation parameter values are 
Am^i = 7Mtofs X IQ-^eV^, taii^ 0^2 = OASltoml 
and sin^ 6*13 = 0.010lg'Q34. The contours are nearly symmet- 
ric about tan^ 6*12 = 1, but the best-fit values for tan^ 6*12 > 1 
are slightly disfavored over those for tan'^ 6*12 < 1, with 
Ax^ =0.8. Assuming CPT invariance, the oscillation pa- 
rameter values from a combined analysis including constraints 



TABLE II: Contributions to the systematic uncertainty in the neutrino 
oscillation parameters Am|i, 612, and 613 for the earlier / later pe- 
riods of measurement, denoted in the text as Period 1 / Period 2 & 3. 
The overall uncertainties are 3.5% / 4.0% for Period 1 / Period 2 & 3. 





Detector-related (%) 


Reactor-related (%) 




Amii 


Energy scale 


1.8/1.8 


Fe-spectra 1321 


0.6 / 0.6 


Rate 


Fiducial volume 


1.8/2.5 


!7e -spectra [241 


1.4/1.4 




Energy scale 


1.1 / 1.3 


Reactor power 


2.1/2.1 




Lcut{Ep) eff. 


0.7/0.8 


Fuel composition 


1.0/1.0 




Cross section 


0.2/0.2 


Long-lived nuclei 


0.3/0.4 




Total 


2.3/3.0 


Total 


2.7/2.8 
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FIG. 6: Prompt energy spectrum of the Fe events in the low-energy 
region for all data taking periods. Bottom panel: data together with 
the best-fit background and geo Ve contributions. The fit incorporates 
all available constraints on the oscillation parameters. The shaded 
background and geo Ve histograms are cumulative. Middle panel: 
observed geo spectrum after subtraction of reactor Fe's and other 
background sources. The dashed and dotted lines show the best-fit 
U and Th spectral contributions, respectively. The blue shaded curve 
shows the expectation from the geological reference model of JlTh . 
Top panel: the energy-dependent selection efficiency. 



from solar neutrino experiments are tan^ 9i2 ~ 0.437!to o26' 
Am|i = 7.53to;ig x IQ-^ cV^ and sin^ 6*13 = 0.023to:oi5. 
A global analysis including also constraints on ^13 from accel- 
erator and short-baseline reactor neutrino experiments yields 



tan2 6'i2 = 0.436 



0.029 
-0.025 
+0.002 
0.002 



Ar 



'21 



7.53 



+0.18 
-0.18 



X 10^5 eV 



and sin^ 6*13 = 0.0231'^ ;;'o2- The fit values for the different 
combinations are summarized in Table HIH Figure |4]shows the 
extracted confidence intervals in the (tan^ 612, Am^i) plane 
with and without the 6*13 constraint. 
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FIG. 7: (a) Confidence level (C.L.) contours for the observed number 
of geo events. The small shaded region represents the prediction 
of the reference model of IitIi . The vertical dashed line represents 
the value of (A^u — NTh)/{Nv + Ntyi) expected for a Th/U mass 
ratio of 3.9 derived from chondritic meteorites, (b) Ax^-profile from 
the fit to the total number of geo Fe events, fixing the Th/U mass ratio 
at 3.9. The grey band represent the geochemical model prediction, 
assuming a 20% uncertainty in the abundance estimates. 



TABLE III: Summary of the fit values for Am|i, t&n^ 612 and 
sin^ 6*13 from three-flavor neutrino oscillation analyses with various 
combinations of experimental data. 



Data combination 


Amil 


tan^ 612 


sin^ 013 


KamLAND 
KamLAND + solar 
KamLAND + solar + 


7 ^4+0-19 
'■04_o.l8 

7 t:q+0.19 
'■O'J-0.18 

5i3 7.53t°:l«8 


n 481+0-092 
U.48i_o 

n 4Q7+0.029 

U.4JD_o.o25 


U.U1U_Q 

n23+° °i^ 

U.UZO_Q Qig 

U.UZO_Q 002 



The KamLAND data illustrates the oscillatory shape of re- 
actor Fe's arising from neutrino oscillation. The ratio of the 
background- and geo-Fe -subtracted reactor 77e spectrum to the 
no-oscillation expectation is shown in Fig.|5]as a function of 
Lq/E, where Lq = 180 km is the flux-weighted average re- 
actor baseline. The improved determination of the geo 17^ flux 
resulting from the addition of the reactor-off data makes the 
second peak at Lq/E = 7Gkm/McV more evident than in 
previous analyses. 
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For the geo flux measurement we incorporate all avail- 
able constraints on the oscillation parameters. The insets in 
Fig.|3]detail the observed spectra in the low-energy region for 
each data taking period. Figure |6] shows the measured geo 
77e event spectrum after subtracting the best-fit reactor and 
background spectra. The best-fit to the unbinned data yields 
1 16 and 8 geo T7p's from U and Th decays, respectively. The 
joint confidence intervals for the sum N\] + A^xh and the asym- 
metry factor {N\j — iVTh)/(-^u + -^Th) are shown in Fig.|7] 
This result agrees with the expectation from the geological 
reference model of UtIi. We obtained an upper limit of <19 
(90% C.L.) in the Th/U mass ratio, indicating the separation 
of U and Th Fg's. Assuming a Th/U mass ratio of 3.9, as pre- 
dicted by the geochemical model of ITTIl from the abundances 
observed in chondritic meteorites, the total number of U and 
Th geo events is 116l;27' ^"^^ a Ax^-profile as shown in 
Fig-Elb). This result corresponds to an (oscillated) Vf. flux of 
3.4lo;g X 10'^ cm^^s"^ at KamLAND, or a total antineutrino 
flux including afl flavors of 6.2l^ g x lO^cm^^s^^. From 
the Ax^-profile (Fig.lTjb)), we find that the null hypothesis is 
disfavored with a p-value of 2 x 10^®. 

The KamLAND data also tests the hypothesis of a natural 
nuclear reactor in the Earth's core ll33ll assuming a constant 
power output over the duration of the experiment. The oscil- 
lation parameters are constrained from the solar, accelerator, 
and reactor neutrino data, while the contributions from geo- 
logical reactor Fe's and from U and Th geo Fg's are allowed 
to vary. The fit gives a limit on the geological reactor power 
of <3. 1 TW at 90% C.L. «3.7 TW at 95% C.L.), an improve- 
ment of a factor of L7 over the previous KamLAND result JH, 
due primarily to the reduction of the commercial reactor i/g 
background in Period 3. 



VII. CONSTRAINTS ON EARTH MODELS 

While the mantle is the most massive layer of the Earth's 
interior, its chemical composition is still uncertain. A quan- 
titative estimate of the heat production by radiogenic compo- 
nents is of particular importance for understanding dynamic 
processes such as mantle convection. Indeed, precisely how 
the mantle convects is still not fully understood, and contro- 
versy remains as to whether two-layer convection or whole- 
volume convection provides a more accurate description. In 
this work, we carry out a comparison of existing Earth models 
using the KamLAND geo T/g data on the basis of simple but 
appropriate assumptions. 

The crustal contribution to the flux at KamLAND can be es- 
timated from compositional data through rock sampling UtI . 
Since current Earth models predict that the lithophiles U and 
Th are absent in the core, for a first approximation of the radio- 
genic heat, we attribute any excess above the crustal contribu- 
tion to U and Th uniformly distributed throughout the mantle. 
Under these generic assumptions, the measured KamLAND 
geo flux translates to a total radiogenic heat production of 
11.2^g'^TW from U and Th. This calculation accounts for 
crustal uncertainties of 17% and 10% for U and Th, respec- 
tively, including correlated errors as suggested in iImI . To 
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FIG. 8: Geo flux versus radiogenic heat from the decay chains of 
U and Th. The measured geo Ve flux (gray band) is com- 
pared with the expectations for the different mantle models from 
cosmochemical (3q1 , geochemical fllll . and geodynamical (stIi es- 
timates (color bands). The sloped band starting at 7TW indicates 
the response to the mantle flux, which varies between the homo- 
geneous and sunken-layer hypotheses (solid lines), discussed in the 
text. The upper and lower dashed lines incorporate the uncertainty in 
the crustal contribution. 



parameterize the planetary-scale energy balance, the fraction 
of the global heat production from radioactive decays, the so- 
called "Urey ratio", is introduced. Allowing for mantle heat 
contributions of 3.0 TW from other isotope decays lfl2l l35l . 
we find that the convective Urey ratio, the contribution to the 
Urey ratio from just the mantle, is between 0.09 and 0.42 at 
68% C.L. This range favors models that allow for a substan- 
tial but not dominant contribution from the Earth's primordial 
heat supply. 

Several established estimates of the BSE composition give 
different geo Ve, flux predictions. Reference 1|38J categorizes 
the models into three groups: geochemical, cosmochemical, 
and geodynamical. Geochemical models ifTHl . such as the 
reference Earth model of iflTll, use primordial compositions 
equal to those found in CI carbonaceous chondrites, but al- 
low for elemental enrichment by differentiation, as deduced 
from terrestrial samples. Cosmochemical models ll36l assume 
a mantle composition similar to that of enstatite chondrites, 
and yield a lower radiogenic abundance. Geodynamical mod- 
els 1331, on the other hand, require higher radiogenic abun- 
dances in order to drive realistic mantle convection. 

In Fig [8] the observed geo flux at KamLAND is 
compared with the expectations from these BSE composi- 
tional models assuming a common estimated crustal contri- 
bution 1 17]. The Fe flux predictions vary within the plotted 
vertical bands due to uncertainties in both the abundances of 
radioactive elements in the mantle as well as their distribu- 
tions. The spread of the slope reflects the difference between 
two extreme radiochemical distributions: the "homogeneous 
hypothesis" in which U and Th are assumed to be distributed 
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uniformly throughout the mantle, and the "sunken-layer hy- 
pothesis", which assumes that all of the U and Th below the 
crust collects at the mantle-core interface. While the statis- 
tical treatment of geological uncertainties is not straightfor- 
ward, assuming Gaussian errors for the crustal contribution 
and for the BSE abundances, we find that the geodynamical 
prediction with the homogeneous hypothesis is disfavored at 
89% C.L. However, due to the limited statistical power of the 
data, all BSE composition models are still consistent within 
-2cr C.L. 



VIII. CONCLUSION 

An updated KamLAND measurement of T^e's was pre- 
sented. This data benefits from the significant reduction of 
reactor Fe's due to the long-term shutdown of commercial nu- 
clear reactors in Japan. The geo Ve_ flux estimate is signifi- 
cantly improved by the reactor-off data. Likewise, the reactor 
neutrino oscillation parameters are also better determined due 
to the reduction of uncertainties in the geo Vf. flux and the 
rates of other backgrounds. Including constraints on 6*13 from 
accelerator and short-baseline reactor neutrino experiments, a 
three-flavor analysis of solar and KamLAND data gives fit val- 
ues for the oscillation parameters of tan^ 612 = 0-436]';Q g25, 
Amii = X 10-5 eV^ andsin^ 6I13 = 0.023tg:^^2- 



Assuming a chondritic Th/U mass ratio of 3.9, we observed 
II6I27 8^° events, which corresponds to a geo 77e flux of 
3.4+o;g X 10^ cm-2g-i KamLAND. The observed rate is 
in agreement with the predictions from existing BSE composi- 
tion models within ^2a C.L. Currently, the ability of discrim- 
inating between models is limited by the experimental un- 
certainty. In the future, improved measurements with higher 
statistics and lower background can be achieved by larger de- 
tectors distant from commercial reactors. Multi-site flux data 
at a combination of crustal and oceanic geological sites would 
be able to estimate the crustal contribution from a statistical 
correlation analysis and constrain mantle abundances more 
stringently. 
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